Measurement and analysis of the pp — y pp~f reaction at 310 MeV 
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The pp — > pp"/ reaction has been studied at a beam energy of 310 MeV by detecting both final 
protons at the PROMICE-WASA facility and identifying the photon through the resulting missing- 
mass peak. The photon angular distribution in the center-of-mass system and those of the proton- 
proton relative momentum with respect to the beam direction and to that of the recoil photon were 
determined reliably up to a final pp excitation energy of E pp ~ 30 MeV. Except for very small E pp 
values, the behavior of these distributions with excitation energy is well reproduced by a new refined 
model of the hard bremsstrahlung process. The model reproduces absolutely the total cross section 
and its energy dependence to within the experimental and theoretical uncertainties. 

PACS numbers: 25.40.Ep; 25.20.-x; 13.60.-r 



I. INTRODUCTION 

The classic motivation for measuring the emission 
of bremsstrahlung in nucleon-nucleon collisions was the 
study of the off-shell behavior of the associated elastic 
scattering amplitude, although it is now known that off- 
shell effects cannot be measured, even in principle. Nev- 
ertheless, the bremsstrahlung reaction can provide a win- 
dow into the underlying dynamical mechanisms that gov- 
ern the NN interaction and the electromagnetic currents 
of nucleons and mesons alike. This is especially true for 
the hard part of the bremsstrahlung spectrum, where the 
photon takes a large fraction of the available center-of- 
mass energy. In this region the A(1232) isobar may also 
play some role and the reaction has then much in com- 
mon with meson production. 

Hard bremsstrahlung has been studied in the radiative 
capture in neutron-proton scattering, np — > dr/, to quite 
high energies through the measurement of the inverse re- 
action of deuteron photodisintegration [l|, Q . The energy 
dependence of the cross section provides direct evidence 
that one of the main driving terms is the excitation of an 
S'-wave AN pair that de-excites through an Ml transi- 
tion into dr/. 

The situation is radically different in proton-proton 
collisions. The analogous Ml transition is forbidden for 
pp —y ppr/ by angular momentum and parity conservation 
when the two protons emerge with very low excitation en- 
ergy Epp = W pp — 2m p , where W pp is the total energy 
of the final pp subsystem in its rest frame Q ■ There are 
therefore significant cancelations among the large ampli- 
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tudes in the pure S'-wave diproton limit, so that any A 
effect must enter in more subtle ways. Furthermore, the 
El transition is generally suppressed by the vanishing of 
an electric dipole operator for the proton pair. It is there- 
fore to be expected that the bremsstrahlung production 
rate should be much lower in pp collisions than in np. 

One technique used to investigate the emission of hard 
bremsstrahlung in proton-proton collisions is the photo- 
disintegration of a pp pair in 3 He. Events where two 
fast protons emerge from 7 3 He — > ppn are interpreted 
in terms of an interaction on a diproton, with the neu- 
tron merely appearing as a spectator 0Q. Such data 
show little evidence for the presence of an intermediate 
AN pahjCertainly much less than for those with fast pn 
pairs [5], Q. However, because the capture on pn spin- 
triplet pairs is so much stronger, the pp data extracted 
in this way may be contaminated by final state interac- 
tions, possibly involving np charge exchange. This can 
only be checked through direct pp —> ppr) measurements. 

Proton-proton bremsstrahlung has been studied in sev- 
eral experiments but, in general, these were undertaken 
by detecting the emerging protons in pairs of small coun- 
ters, often placed on either side of the beam direc- 
tion @, [1] , which has led to the low E pp region being es- 
pecially poorly sampled. The geometric acceptance was 
much increased in a series of refined KVI experiments at 
190 MeV [ildtl, but even here the low E pp region was not 
favored. Whereas the COSY time-of-flight spectrometer 
also has wider acceptance, the data obtained at 293 MeV 
have only limited statistics and no attempt was made to 
evaluate the cross section as a function of E pp [111 ]. 

Data on the hardest part of the pp bremsstrahlung 
spectrum were also obtained at the COSY-ANKE mag- 
netic spectrometer by selecting the two final protons with 
Epp < 3 MeV [H, [H| ■ A proton beam energy range from 
353 to 800 MeV was investigated but only for CM pho- 
ton angles 6 1 , where cos# 7 > 0.95. The results reveal 
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a broad peak in the cross section at an energy around 
650 MeV with a FWHM 220 MeV. This suggests the 
possible influence of intermediate AN pairs, though not 
necessarily in a relative S-wave. 

Much higher statistics were obtained over a wider 
range of pp excitation energies and photon angles at the 
PROMICE-WASA facility at Uppsala. The experiment 
was carried out at a single beam energy of T p = 310 MeV 
and results were recently published for E pp < 3 MeV [3] ■ 
These data are completely consistent with those from 
ANKE at 353 MeV over the small-angle domain cov- 
ered by the ANKE experiment [HI, [l3j]. However, it is 
clear from this comparison that a reliable decomposi- 
tion into multipoles requires data over a wide angular 
range. The low E pp data from Uppsala were interpreted 
as indicating the dominance of the El and M2 multi- 
poles with no evidence for any important E2 contri- 
bution, in contrast to theoretical expectations [l5|. The 
purpose of the present paper is to extend the analysis up 
to E pp w 30 MeV in order to test theoretical models over 
a wider range of excess energies. 

A state-of-the-art model has recently been developed 
that for the first time describes successfully proton- 
proton bremsstrahlung in the hundred MeV range 
[111 . This model, which is hereinafter denoted as HN, 
is summarized in Sec. [TXJ In this approach the photon 
is coupled everywhere to a relativistic pp scattering am- 
plitude in a way that ensures consistency with gauge in- 
variancc. Although this reproduces very well the detailed 
KVI pp —} ppj measurements at 190 MeV @, it is pos- 
sible that at 310 MeV the tail of the A might have some 
influence. In this context it should be noted that the 
minimal inclusion of the A isobar [l8[ (see also Ref . [lH ) 
improves the theoretical description of the 280 MeV TRI- 
UMF data ||. 

The experimental approach used in this work is identi- 
cal to that employed at PROMICE-WASA for pion pro- 
duction [IjJ |2fJ and so Sec. IIIII and Sec. IIVI merely pro- 
vide outlines of the salient points of the method and the 
data analysis, respectively. The results given in Sec. fVl 
show that, away from the region of small E pp values, 
where there can be significant cancelations between dif- 
ferent contributions, the theory of Sec. [TT] works remark- 
ably well. It describes the photon angular distribution 
and those of the diproton relative momentum in differ- 
ent Epp intervals as well as the energy dependence of the 
pp — > ppj total cross section. The fact that the the- 
ory reproduces the absolute normalization of these high- 
momentum-transfer data to within the experimental and 
theoretical uncertainties is striking. However, the the- 
oretical predictions of the photon angular distributions 
obtained without intermediate AN contributions are bet- 
ter at low Epp than those that include them. This brings 
into question whether the present simplified treatment 
of these isobar contributions is acceptable. Our conclu- 
sions and suggestions for further work are to be found in 
Sec. EH 



II. THEORETICAL MODEL 

The experimental data presented in this paper repre- 
sent the most complete measurement of bremsstrahlung 
in proton-proton collisions at an energy so far above 
threshold. We therefore compare them with a re- 
fined mod el [lq that has recently been successfully ap- 
plied [H Qjto describe both the TRIUMF data at 
280 MeV [1] and the high-precision KVI data [lj at the 
lower proton energy of T p = 190 MeV. This solved a long- 
standing discrepancy between experiment and the, then, 
existing theory. The novel approach is derived within 
a quantum field-theory formalism by coupling the pho- 
ton everywhere possible to an underlying two-nucleon T- 
matrix that is derived from a relativistic NN scattering 
equation. The basic idea of the method is that intro- 
duced by Haberzettl, Nakayama, and Krewald [2l[ for 
pion photoproduction, based on the field-theoretical ap- 
proach of Haberzettl (22|. The model accounts for the 
important interaction current in the NN bremsstrahlung 
reaction in a manner that is consistent with the gen- 
eralized Ward-Takahashi identity (WTI), which ensures 
gauge invariance at the microscopic level. This feature is 
absent from all earlier models. 

Following Ref. one starts from the nucleon- 

nucleon T-matrix determined by the relativistic Bethe- 
Salpeter (BS) equation 

T = V + VG T = V + TG V , (1) 

where V represents the driving two-nucleon potential. 
The two-nucleon propagator, Go = S1S2, describes the 
intermediate propagation of two free non-interacting nu- 
cleons (with individual Feynman propagators Si, i = 1,2) 
sharing the given fixed reaction energy. This relativis- 
tic four-dimensional equation is then reduced in a co- 
variant manner to the three-dimensional Blankenbecler- 
Sugar (BbS) equation [23|, |2J] by replacing the propaga- 
tor Go by Go — > go where go restricts the intermediate 
two nucleons to be on their mass-shells in a manner that 
preserves the (relativistic) unitarity of the equation. 

The driving potential V used here is based on the one- 
boson-exchange model developed by the Bonn group [24| , 
which contains nucleonic and mesonic degrees of freedom. 
In addition to reproducing the low-energy pp scattering 
data and the deuteron properties, the resulting NN in- 
teraction fits the NN phase-shifts up to the threshold 
for pion production. This version is used, rather than 
a more modern potential, because the necessary interac- 
tion current that is fully consistent with this potential is 
already available from the work of Ref. pjj > where it was 
shown to be crucial in resolving longstanding theoretical 
issues with the KVI data 

By coupling the photon to the system of two interact- 
ing nucleons, it can be shown, again following Ref. [lj|, 
that the resulting bremsstrahlung amplitude may be 
written as 

M» = {Tg + l)J»(l+g T) , (2) 
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FIG. 1. (a) Basic photon production current J M used to describe the pp — > ppj reaction in the present work; analogous photon 
couplings along the lower nucleon line are omitted. N denotes an intermediate nucleon, and M incorporates all the meson 
exchanges. As indicated by the symbols below the diagrams, the first two describe the nucleonic current, while the meson- 
exchange current is depicted by the third. The fourth diagram contains the NM — > N'y four-point contact current which, 
together with the meson-exchange current, constitutes the interaction current of Eq. ([3]). Diagram (b) shows the photon 
coupling to both intermediate nucleons as subsumed in the dual current d M . 



where the final-state (FSI) and initial-state interactions 
(ISI) are included through the NN T-matrices on the left 
and right, respectively. 

The basic photon production current from the two nu- 
cleons 



= d"G V + VG d" + 



(3) 



contains nucleonic and mesonic terms as well as a four- 
point contact-type term, as illustrated in Fig. [lja). 

The two disconnected nucleonic terms, shown in 
Fig. [TJb) and subsumed in the dual current cZ M , are given 
by 



(4) 



Here is the NNj vertex for nucleon i (= 1,2), Si de- 
notes the propagator of the nucleon i, and 8t represents 
an implied <5-function that ensures that the incoming and 
outgoing momenta of the intermediate nucleon i are iden- 
tical. 

We mention that the dynamical structure of this for- 
mulation takes care of the fact that the translation of the 
three-dimensional BbS reduction to the bremsstrahlung 
reaction must be implemented such that a physical pho- 
ton cannot couple to a nucleon that is on-shcll before and 
after the coupling takes place. For more details regarding 
this non-trivial issue, see Ref. (l6| . 

The of Eq. describes the photon coupling to 
the internal mechanisms of the interaction V, i.e., it cor- 
responds to the interaction current. For a one-boson- 
exchange model of the NN interaction, such as that 
employed here, V^- consists of mesonic and four-point 
contact currents. Unlike the case of proton-neutron 
bremsstrahlung, where there is a large mesonic current 
contribution [25|, [26| , this is to a large extent suppressed 
for proton-proton bremsstrahlung because only neutral 
mesons can then be exchanged. The dominant mesonic 
current contributions that we include arise from the 
anomalous vttj couplings (v — p,u>). These transitions 
are transverse and thus cannot be obtained by simply 
coupling the photon to the underlying NN T-matrix; 
they must be inserted by hand into J M . 

The four-point contact current appears as a conse- 
quence of imposing gauge invariance in the form of the 
generalized WTI on the resulting amplitude. Note that 



the vnj meson-exchange currents have no influence on 
this because they are purely transverse. In general, 
contact-type currents have very complicated microscopic 
dynamical structures [22| that cannot be taken into ac- 
count explicitly at present. Instead, one must revert 
to employing generalized phenomenological contact cur- 
rents, constructed such that the full reaction amplitude 
satisfies the generalized WTI, which is necessary to en- 
sure full gauge invariance at the microscopic level. As 
a consequence, no unique determination of the reaction 
amplitude is possible since the WTI does not constrain 
the transverse part of the amplitude. In the present case, 
the dynamics of the hadron interactions is described in 
terms of phenomenological form factors. The resulting 
phenomenological four-point interaction currents that 
describe the interaction of the photon with this hadronic 
three-point function, therefore, are constructed purely in 
terms of these hadron form factors [H, [TtJ ■ 

In this paper we use our dynamical model in the anal- 
ysis of the pp — > ppj reaction data at a proton inci- 
dent energy of 310 MeV. Although this is well below the 
maximum of A production at about 650 MeV, earlier 
analyses Q3 (see also Ref. [Hj]) of the TRIUMF data 
at 280 MeV [8| show that its inclusion can improve the 
agreement with data in certain geometries. We therefore 
investigate the effect of introducing the A in a minimal 
fashion, following the application section of Ref. 16| , by 
implementing the A contributions in J M at the tree-level. 
For this purpose two more terms, analogous to the first 
two on the r.h.s. of Fig. [Ha), are added, with the A 
resonance replacing the intermediate nucleon N. This 
A resonance current has no bearing on gauge invariance 
because the A./V7 transition vertex is purely transverse. 
However, in a full NN ^ AiV coupled-channels ap- 
proach, in addition to the tree-level A resonance current 
considered here, there will also be additional box-type 
contributions with intermediate AN and AA pairs that 
produce purely transverse five-point contact-type contri- 
butions to the interaction current [l6| . At this stage, 
therefore, the present minimal tree-level inclusion of the 
A currents should be considered only exploratory. 

The Bonn potential employed in the present study for 
generating the nucleon-nucleon T-matrix is given in mo- 
mentum space; it is therefore non-trivial to include the 
Coulomb interaction. Coulomb effects have been investi- 



4 



gated inpp bremsstrahlung in the past (27j . However, the 
associated distortions are mainly relevant at very small 
pp invariant masses, a regime which has not been well 
sampled in most of the earlier experiments. 

In order to test the influence of the Coulomb inter- 
action over the wider acceptance of the present exper- 
iment, we also consider the NN interaction based on 
the Paris potential [28]. The Paris work was carried 
out in coordinate space and includes fully the Coulomb 
interaction (2?j but only within the framework of the 
non-relativistic Lippman-Schwinger equation. We have 
therefore formally transformed this into the relativistic 
BbS equation by a proper redefinition of the potential, 
through the so-called minimal relativity factor [29[ , in or- 
der to be able to use this interaction consistently within 
the present relativistic approach. The transformed inter- 
action reproduces the same nucleon-nucleon observables 
for relativistic kinematics as the original one for non- 
relativistic kinematics. 

One shortcoming in the present approach for incorpo- 
rating the Paris potential is that, for simplicity, we have 
retained the production current J M calculated from the 
Bonn potential. As a result, the consistency of the initial 
and final state interactions with the production current 
J M is lost but, for the purpose of checking the Coulomb 
effects, this inconsistency is not of major concern. 
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FIG. 2. Distribution in the square of the missing mass in 
the pp —¥ ppX reaction presented in units of the neutral pion 
mass. Clear peaks are seen, arising from the pp — > ppw° and 
pp — > pp'y reactions, sitting on a slowly varying background. 



III. EXPERIMENT 

The data of the present experiment were obtained 
at the The Svedberg Laboratory in Uppsala, where a 
48 MeV proton beam from the cyclotron was injected 
into the CELSIUS ring ||cj], accelerated to 310 MeV and 
then stored. An average beam-on intensity of 3 mA was 
achieved during an experimental data-taking period of 
approximately 100 hours. 

The measurements were carried out at the PROMICE- 
WASA facility !3J| and results at small proton-proton ex- 
citation energies have already been published [l4j ]. Fur- 
thermore, the pp —¥ ppj data were obtained simultane- 
ously with those on pp — > ppir 19|, whose results were 
reported in greater detail in Ref. [20( . Since the detector 
assembly and the measurement techniques were identical 
in the two experiments, and the experimental procedures 
and data analysis differed only in minor details, the de- 
scription will here be kept quite brief. 

An internal gas-jet hydrogen target, with a density of 
about 2 x 10 14 cm~ 2 , was used in conjunction with the 
stored proton beam. By operating an electron cooler 
throughout the experiment, the background was reduced 
significantly and the counting rate increased through an 
improvement in the beam-target overlap. 

Even though the PROMICE-WASA facility was 
equipped to detect high energy photons, in the 
bremsstrahlung study reported here, only protons were 
measured in the final state. After exiting the scattering 
chamber, the protons passed through a forward window 



counter (FWC), a tracker, a forward trigger hodoscope 
(FTH) and usually stopped in a forward range hodoscope 
(FRH). The four-quadrant scintillator of the FWC elim- 
inated most of the beam-halo background. In order to 
be accepted by the main trigger, coincident protons must 
appear in different quadrants. Events with protons in the 
same quadrant were allowed by a secondary trigger but, 
in accord with Monte Carlo expectations, these were very 
few in number and were not considered in the subsequent 
analysis. 

Information on the proton angles was extracted from 
the FTH and, even more precisely, from the tracker. 
Events with polar angles between about 3° and 22° were 
recorded. As described fully in Ref. [2(|, the energy as- 
sociated with a proton track was deduced from a combi- 
nation of the calculated angle-dependent range up to the 
entrance of the stopping scintillator and the measured 
light output of that detector. A few protons stopped 
in one of the thin dead regions between the scintillator 
planes and these were then assigned the energy corre- 
sponding to the midpoint of the dead layer. 

As an extra check on the particle identification, it was 
further required that both protons of an accepted event 
penetrate at least into the second layer of the FTH, which 
consists of 24 spiral scintillator segments. The minimum 
energy of each proton was therefore 38 MeV. This con- 
dition meant that all coincident pairs of protons stopped 
in the second FRH scintillator or earlier so that there 
was effectively no high-energy limitation imposed by the 
design of the apparatus. 
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FIG. 3. (Color online) Angular distribution of the final pro- 
tons from the pp —¥ pp^f reaction in the laboratory (open 
circles) compared to Monte Carlo simulations based upon 
the HN model [l£|, using the Bonn potential as input (red 
solid curve). Similar predictions obtained with a phase-space 
model (blue dashed curve) are also shown. Both sets of pre- 
dictions were normalized to the total number of events. 



In the missing-mass distribution of the pp — > ppX° 
reaction shown in Fig. [2J there are two clear peaks cor- 
responding to X° — 7 and X° — ir°, with very little 
overlap. Before making any detailed cuts, these peaks 
contained in total 66,521 pp'y and 861,449 ppir candi- 
dates. The exclusion of events affected by the detector 
gaps, and those where the proton time difference fell out- 
side a 65 ns band, eliminated 7.3% and 1.5% of these, 
respectively. There is only a small (« 5%) background 
under the 7 peak that arises mainly from the rescattering 
of one of the protons from a pion-production reaction. 
The maximum polar angle of protons from ir° produc- 
tion depends sensitively upon the proton beam energy 
T p . A measurement of this angle, which was close to 18°, 
showed that T p = 309.7 ± 0.3 MeV. 

The width of the 7 peak is a(M%) = 0.056M2 f) . By 
retaining only events at a little over the two FWHM level, 
namely \M\ /M% 1 < 0.137, to a good approximation this 
cut compensates for the neglect of the small background 
contribution [14j . 

The angular distribution of the protons in the labora- 
tory system for the selected pp — > ppr/ events is shown in 
Fig. [3] In spite of a slight but significant misalignment 
between the beam and the detector axes, the angular 
cutoffs at both small and large angles are quite sharp 
and very well reproduced by the Monte Carlo simulation 
that used the Bonn pp potential in the model described 
in Sec.|ll] The phase-space simulation gives a marginally 



FIG. 4. (Color online) Time difference between the two pro- 
tons emerging from the target. The peaks have similar widths 
for photons in the forward (blue triangles) and backward (red 
circles) CM hemisphere, as well as their sum (black crosses). 
Events in the interval ±1.8 ns were retained for the subse- 
quent analysis. 



poorer representation, especially at large angles. In both 
cases the predictions have been normalized to the total 
number of events. 

For each of the emerging protons, a timing signal was 
extracted from the first of the FTH spiral detectors. The 
time at the target position was then estimated using the 
information on the particle energy, the hit position in 
the scintillators, and the time-of- flight. The calibration, 
which was improved over that used in Ref. [2(|, led to 
a distribution for the time difference between the two 
protons with a peak width of 1.1 ns FWHM. This was 
essentially the same for both the forward- and backward- 
going photons, though a correction was introduced to 
compensate for a slight offset of 0.2 ns in the forward case. 
Cuts at ±1.8 ns applied to the data of Fig. [4] reduced 
the number of accidental coincidences to less than 1% so 
that it was then justified to employ a kinematic fitting. 
This was achieved by adjusting the energies of the two 
protons to give zero missing mass. Taken together with 
the sharper time difference cut, this reduced the number 
of candidates to the 58,026 pp^f events that were used in 
the subsequent analysis. 



IV. ANALYSIS 

In the case of a production reaction like pp —> pp'y, the 
unpolarized cross section is a function of four indepen- 
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dent variables. The standard set chosen for the analysis 
consists of 

• E pp : the excitation energy in the final pp system, 

• 6* 7 : the CM production angle of the photon, 

• 9 q : the CM polar angle of the pp relative momen- 
tum q with respect to the beam direction, 

• (f q : the azimuthal angle between q and the photon 
momentum. 

Other variables, such as the laboratory proton angle that 
was used in the construction of Fig. [31 can be expressed 
in terms of these four quantities. 

In order to convert the observed numbers of events into 
cross sections, knowledge of the detector acceptance is 
needed in the four-dimensional space. This was achieved 
using a Monte Carlo simulation, where the detector sys- 
tem was described in great geometric detail. Identical 
cuts were then placed on the simulated and experimen- 
tal events. Frequently only phase-space was used in the 
simulations but, in principle, the acceptances might de- 
pend significantly on the actual reaction probability. Es- 
timates were therefore made, not only for simple phase 
space, but also for a realistic reaction matrix, assumed to 
be represented by the HN model 16]. In the latter case 
the program interpolated within a lookup table of the re- 
action matrix yielded by this model in the four standard 
variables. Values of the acceptance in the {E pp , 7 ) space 
are shown in Fig. 

It should first be noted that the PROMICE-WASA 
detector only registered protons with laboratory angles 
less than 22° and at 310 MeV this means that only the 
region E pp < 42 MeV was sampled. Although at low 
E pp the acceptance could be quite large, being up to 
60%, this decreased to much lower values at higher E pp 
and small 7 . This is illustrated in Fig. [5] by showing 
the acceptance as a function of E pp in four ranges of 
the photon (CM) angle. These are, respectively, the 
backwards —1 < cos6> 7 < —0.8 (bw), the backwards 
central —0.8 < cos6> 7 < (bwc), the forwards central 
< cos# 7 < 0.8 (fwc), and the forwards 0.8 < cos$ 7 < 1 
(fw) regions. For photons emitted in the forward hemi- 
sphere, the recoiling protons are slower and a greater 
fraction emerge at larger angles than allowed for in de- 
sign of the PROMICE-WASA detector, and this leads to 
a more severe cut at high E pp . Protons from these events 
are also more likely to be distorted by secondary interac- 
tions. On the other hand, it also means that the beam- 
pipe effect kicks in at lower E pp , which is also clearly seen 
in Fig. [SJ 

Due to the identical nature of the protons in the en- 
trance channel, the pp — > ppj cross section is symmetric 
in the CM system around 6* 7 = 90°. The effects of the 
variation of the acceptance with d 1 at large E pp can also 
be seen in Fig. [51 which shows the cross sections extracted 
as functions of E pp in the same four regions of cos # 7 used 
in Fig. [5] In all cases the data were terminated when the 
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FIG. 5. (Color online) Monte Carlo estimates of the percent- 
age acceptance for pp — > pp'y at T v = 310 MeV obtained using 
the HN model pi| with the Bonn potential. The results are 
divided into four regions in the photon CM angle # 7 defined in 
the text, namely the backwards (red circles), the backwards 
central (blue triangles), the forwards central (magenta stars), 
and the forwards (green inverted triangles). 



estimated acceptance dropped below 2%, and this oc- 
curred much earlier for small values of B 1 . Except at the 
edges of the acceptance, the cross sections deduced using 
a phase space model to evaluate the acceptance differed 
only marginally from those obtained on the basis of the 
dynamical model. 

The crucial forward/backward symmetry is clearly re- 
spected to within the uncertainties for E pp < 20 MeV but 
between 22 and 26 MeV there is some deviation, which 
is more apparent in the angular distributions to be pre- 
sented in Sec. [VI On general grounds one would expect 
the data from the forward photon hemisphere to be less 
reliable because the associated protons are less energetic 
and can emerge at larger angles. The statistics in the 
backward hemisphere are also much larger. 

The integrated luminosity of 340 nb _1 was derived 
from a comparison of elastic proton-proton scattering re- 
sults measured in parallel with tabulated cross sections, 
as described in Ref. [2(| ■ Due to the large pre-scaling 
factor associated with the pp trigger used, an error bar 
of about 10% must be associated with this value. This 
includes also effects connected with the evaluation of the 
proton acceptance in the apparatus and any uncertainty 
in the pp database used in the comparison. 

Of the other systematic uncertainties discussed in 
Ref. [20( | . proton rescattering in the detector material 
might contribute 2%, as might the treatment of the back- 
ground under the 7 peak. Although the PROMICE- 
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FIG. 6. (Color online) Cross section for pp — ► pp*y at 
T p = 310 MeV in arbitrary units, obtained using the HN 
model [l6j with the Bonn potential to evaluate the accep- 
tance. The experimental data are divided into the same four 
angular regions as in Fig. [5] Data were only plotted when the 
acceptance was estimated to be above 2%. 



WASA geometric acceptance is very good, the extrap- 
olation to unexplored regions and its dependence upon 
reaction models can give up to 3%, though this depends 
upon the value of E pp . The known systematic uncer- 
tainty is therefore judged to be « 15% overall. How- 
ever, despite the care taken with the calibrations and 
the evaluation of the acceptance, the forward/backward 
symmetry is not completely respected at high E pp , as ev- 
idenced by the divergence between the fwc and bwc data 
in Fig. [6] We therefore cannot exclude larger systematic 
uncertainties even in the backward photon hemisphere 
for E pp > 30 MeV. 



V. RESULTS 

Over 58,000 kinematically well-defined pp — > ppry 
events are available for analysis in terms of the four- 
dimensional differential cross section, as described in 
Sec. IIVI In the present paper only one single differential 
and three double differential distributions are presented. 
Data points are shown if the acceptance at this point is 
estimated to be larger than 2%. Only the statistical un- 
certainties are shown explicitly by error bars and these 
do not include the ~ 15% overall systematic effects. The 
azimuthal dependence of the data can be quite strong 
but, since this seems mainly to be a reflection of the ac- 
ceptance, it is not further investigated. The data would 



allow explorations of higher dimensionality, but further 
guidance from theory would be necessary to exploit this 
fruitfully. 

The center-of-mass differential cross section in the 
photon angle is presented in Fig. [7] averaged over 3 MeV 
bins in the pp excitation energy from 0-3 MeV to 39- 
42 MeV, with the upper end of each interval being in- 
dicated in the relevant panel. The angular cuts on the 
data clearly reflect the acceptance dependence presented 
in Fig. [SJ 

As already reported [l4j . the data for E pp < 3 MeV 
show a strong minimum at d 1 — 90° and an almost pure 



cos" 



9 7 behavior. The level of this minimum rises as E pp 



/ 7 term is gen- 



mcreases. The data show that the cos 4 9 1 
erally small and its strength cannot be determined with 
precision. 

The full dynamical model of Haberzettl and 
Nakayama [l6[ of Sec. [XT] has been evaluated using the 
Bonn potential [24| without the A contribution. The 
curves are consistent with the shapes of the angular dis- 
tributions as measured in the backward hemisphere for 
all except the lowest E pp bin. The strengths are also 
well described, especially in view of the uncertainties in 
the absolute scales of both the theory and experiment. 
A minimum is predicted at 90° for all energy bins but 
for Epp < 3 MeV this is not sufficiently deep and there 
seems to be no sign there of the leveling off near the for- 
ward/backward directions expected from the theoretical 
model. 

The theoretical predictions are, of course, sensitive to 
the assumptions in the model. Thus, when the "mini- 
mal" inclusion of the A contribution is switched on in 
the calculation, the effects are surprisingly large and the 
agreement with the data is much poorer, especially at 
low Epp where the central minimum is largely absent. 
The shapes are far less changed at high E pp and the data 
there can be well reproduced if the predictions are scaled 
by a factor of « 0.8. 

The description of the angular distributions with the 
Paris potential [2j§] is very similar to that obtained with 
the Bonn potential. The Coulomb effects that are in- 
cluded here are only significant for very low E pp but this 
is also the region where the theoretical model is least 
satisfactory. 

The second angular distribution to be discussed is that 
of the pp relative momentum vector q with respect to 
that of the proton beam in the overall CM frame. This 
is shown in Fig. [8] in the same 3 MeV bins that were 
used in Fig. [7] It is difficult to measure the angles of 
the vector q when its magnitude is small so that any 
apparent deviation from isotropy for E pp < 3 MeV may 
not be significant. 

The data for E pp > 3 MeV show clear evidence of a 



forward dip and the predictions of the HN model [16| on 
the basis of the Bonn potential follow this trend quanti- 
tatively. As can be seen from Fig. [51 in this case there 
is very little difference in the shape of the predictions 
whether the A is included or not. The lack of data for 
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FIG. 7. (Color online) Differential cross section for pp — > pp*y at 310 MeV. The data have been averaged over 3 MeV bins in 
E pp and the upper end of each interval is indicated by the number at the bottom right of the corresponding panel. The curves 
correspond to the absolute predictions of the Haberzettl and Nakayama model 116] of Sec. [TTJ using the Bonn potential [24] . The 
solid (blue) one represents the model without any A contribution whereas for the dashed (red) one the A has been switched 
on. The latter has been included in a minimal fashion, as explained in Sec. [II] 



small | cos 9 q | at high E pp is a consequence of the limited 
acceptance in these regions. 

The corresponding distribution of q with respect to 
the photon direction is shown in Fig. [9] Although this 
is evaluated in the frame of the recoiling pp pair (the 
helicity distribution), little change would be seen if this 
were replaced by the overall center-of-mass frame. Both 
the data and the models display fairly flat shapes, with 
the possible exception of the very low energy bins, where 
the drop in the data for small helicity angle Oh may reflect 
the difficulty in measuring two protons when they emerge 
with similar angles and momenta. 

Note that the data in the three angular distributions 
are different representations of the same 58,026 events, 
so that the cross sections integrated over angle must be 
identical. The "holes" seen at various places of phase 
space in the diagrams indicate possible sources of system- 
atic errors. As a consequence one must conclude that the 
integrated pp —¥ pjry cross section can only be safely ex- 
tracted when E pp < 30 MeV. The energy variation of this 



cross section is shown in Fig. [TU]as a function of E pp . For 
one of the sets of points, only results from the backward 
photon hemisphere (bw + bwc) are used. The other set 
uses in addition the fwc data and the forward/backward 
symmetry is only invoked to derive the data in the very 
forward region. 

The reasonable description of the angular distributions 
by the theoretical model is translated into one of the in- 
tegrated cross section, where the predictions have been 
smeared over the 0.5 MeV bins used in the data presenta- 
tion. In order to ensure agreement with the data above 
10 MeV, the theoretical results that included only the 
nucleonic and meson-exchange current terms were scaled 
by factors of 1.30 and 1.45 for the Bonn and Paris eval- 
uations, respectively. When the AN intermediate states 
are included in the minimal way described in Sec. HH the 
corresponding factors are 0.80 and 0.90, respectively. 

The drawback of using the Bonn potential is immedi- 
ately apparent. The unrealistic spike at very low E pp is 
significantly softened when this is replaced by the Paris 
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FIG. 8. (Color online) The distribution in the angle between the pp relative momentum vector q and that of the incident beam 
direction in the overall CM frame. It should be noted that there may be significant systematic measurement uncertainties at 
low E pp . Since the final protons are identical, data are only shown in one hemisphere and the total cross section is obtained 
by summing over this region. The curves correspond to the predictions of the Haberzettl and Nakayama model [16j of Sec. [II] 
obtained using the Bonn potential [24[. The notation is the same as in Fig. [7) 



potential, which includes the Coulomb repulsion in the 
pp system. Otherwise there is little difference between 
the predictions based upon the two potentials. In either 
case the model seems to overestimate the cross section for 
the production of the ^Sq state of the two final protons. 
It must again be stressed that in this region there can 
be delicate cancelations amongst the contributions [3[ . If 
the ^So prediction were reduced, the energy dependence 
might be reproduced, though with a slightly too low over- 
all normalization. 

The inclusion of the effects of the A in an approxi- 
mate way gives a rather similar description of the data 
in Fig. [lOj However, as mentioned already, the A ef- 
fects beyond the tree-level that would arise in a consis- 
tent NN ^ AiV coupled channels approach [l^|, have 
been ignored in the present calculations. Using the Paris 



rather than the Bonn potential induces changes analo- 
gous to those seen in the non-A scenario. 

Phase space does not provide an acceptable descrip- 
tion of the E pp dependence of the integrated cross sec- 
tion shown in Fig. 1101 This is by no means unexpected 
because the spin-parity constraints associated with the 
So final state are not built into such a naive approach. 

The total pp — > ppj cross section integrated up to 
E pp = 30 MeV is er(30) = (0.59 ± 0.09) fib, where the 
statistical error is negligible and the quoted uncertainty 
is purely systematic. In the COSY-TOF measurement 
at 293 MeV [llj there were very few events with E pp > 
60 MeV and a total cross section of a = (3.5±0.3±0.7) fib 
was obtained by extrapolating to the kinematic limit on 
the basis of a phase-space variation. If, instead, the to- 
tal cross section below 30 MeV is estimated, a value of 



10 



20 

10 


40 

w 20 



•3 ° 

w 40 
S & 20 

S o 

£ 40 



T3 20 




20 



3 


6 




1 1 1 


nm..,,' 

i 

i 1 i 


9 


12 


15 


1 

1 1 1 


i 

i i i 


i 

i l i 


18 


21 

■,•••••■••••«■ 


24 

■ ; : i ■ ' ' ' 


1 

1 1 1 


1 

l 1 l 


i 

i 1 i 




I.l^T"l' 




27 

1 1 L "i— 1 


1 1 1 — ■ 


33 

i 

i I i , 




■ " 


kt 


36 

i 1 i 


* k* 

- i 

1 39 

i I i 


i I i 



0.5 1.0 0.5 1.0 

cos8 h (deg) 



0.5 



1.0 



FIG. 9. (Color online) The distribution in the helicity angle between the pp relative momentum vector q and that of the photon 
in the frame of the recoiling pp pair. The notation is the same as in Fig. [7] 



ct(30) = (0.54 ± 0.12) is found. The agreement be- 
tween the COSY-TOF result and ours lies well within 
the error bars. 

In order to illustrate the variation of the angular de- 
pendence of the photon with excitation energy, the dif- 
ferential cross section data of Fig. [7] have been fitted by 



^^ = 2^ a2 " P2 " (COS ^ ) - (5) 

7 n=0 

Figure QT] shows the ratio 02 /ao as a function of E pp 
compared to the predictions of the present calculation 
based on the HN model [l6[ . Neither the prediction with 
or without the A contribution can describe the rise in 
0-2 1 o-q at low Epp and the inclusion of the A is particularly 
disappointing in this region. 



VI. CONCLUSIONS 

We have presented here detailed measurements of 
bremsstrahlung production in proton-proton collisions at 
a beam energy of 310 MeV. The differences between 
the data extracted from the forward and backward pho- 
ton hemispheres increases for excitation energies above 
about 20 MeV and, in such cases, the backward data 
are more reliable because of the faster protons and the 
much higher acceptance. Using the forward/backward 
symmetry of the reaction, full acceptance was achieved 
up to an Epp w 30 MeV and some information obtained 
even close to the kinematic limit of 42 MeV imposed by 
design of the apparatus. The big advantage of this ex- 
periment compared to others undertaken above the pion- 
production threshold is the large acceptance coupled with 
good statistics. It is therefore not surprising that the 
values extracted for the cross sections depend very little 
whether one uses phase space or the HN model to esti- 
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FIG. 10. (Color online) Cross section for the pp — >• ppr/ 
reaction measured at T p — 310 MeV integrated over all 
angles. The (black) circles were obtained by taking data 
from the backward photon hemisphere and invoking the for- 
ward/backward symmetry. For the (magenta) stars, data in 
the interval —1 < cos# 7 < 0.8 were used in the determi- 
nation. Overall systematic errors of ±15% are not shown. 
The solid (blue) curve represents the evaluation of the model 
described in Sec. [IT] using the Bonn AW input without the 
inclusion of any A contribution, whereas the red dots rep- 
resent an attempt to include this in a minimalist way. The 
(orange) dashed and green chain curves are the analogous es- 
timates based upon using the Paris potential that includes 
the Coulomb repulsion. The curves are scaled by the factors 
shown in the figure. The (maroon) dot-dot-dashed curve il- 
lustrates the energy dependence to be expected from a pure 
phase-space model. It is arbitrarily normalized at around 
5 MeV, i.e., just above the region whose energy dependence 
is largely driven by the strong pp FSI. 
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FIG. 11. (Color online) Ratio of the Legendre coefficients de- 
fined by Eq. (0. Values extracted from experimental data are 
compared to the predictions of the Haberzettl and Nakayama 
model [lfj with and without the minimalist A inclusion (red 
dotted and solid blue lines, respectively). 



mate the acceptance. However, to go to higher E pp one 
would need larger counters than those provided by the 



PROMICE-WASA setup. 

Away from the very small E pp region, the dynami- 
cal model of Haberzettl and Nakayama [16| , whose main 
points are summarized in Sec. UTJ is rather successful in 
describing all the experimental results as functions of E pp 
provided an overall scaling factor close to unity is applied. 
As well as the integrated cross section, these include the 
angular distributions of the photon and those of the pp 
relative momentum with respect to the beam direction 
and to that of the recoil photon. It is very gratifying 
to note that any rescaling of the predictions required 
to achieve this success is well within the combined ex- 
perimental and theoretical uncertainties. The latter are 
clearly very hard to quantify but they must at least en- 
compass the differences between the inclusion or not of 
the A contributions. 

The situation at low E pp is more uncertain because 
there are significant cancelations amongst the driving 
terms [H and the theoretical results are therefore much 
more sensitive to small contributions. Using the Paris 
rather than the Bonn potential in the evaluation of the 
model allows the Coulomb interaction to be included and 
this does smooth the predictions slightly at low E pp . 
However, it must be noted that the switch from Bonn 
to Paris was not done fully consistently. 

The good agreement between theory and experiment 
was obtained without considering any effects that might 
arise from the virtual excitation of the A isobar. Al- 
though the basic model was tuned to describe the 
190 MeV KVI data :9j, by 310 MeV the influence of the A 
might start to be felt. In this context it should be noted 
that the introduction of the A isobar [lj| [l|| seems to 
improve the theoretical description of the 280 MeV TRT- 
UMF data Q. However, we find that the inclusion of the 
A effects in the minimal way described in Sec. |H] actu- 
ally makes the agreement with the shapes of the photon 
angular distribution worse for E pp < 12 MeV and this 
difference is somewhat puzzling. We view it more as an 
indication that A effects are not very well understood 
and that minimal inclusion is not warranted in these in- 
stances, rather than as a measure of theoretical errors. 
Further theoretical work is clearly needed and the results 
might be improved by introducing a phenomenological 
five-point contact current to take account of the AN box 
diagrams, as explained in Ref. (l(ij |. 

At higher E pp the differences are less important and 
the four versions of the model give very similar integrated 
cross sections in Fig. [TUl provided that they are scaled by 
factors that are all fairly close to unity. With the scaling 
factors as shown, the predictions start to differ from the 
more reliable backward-angle data above 20 MeV. How- 
ever, in view of the 15% overall systematic error in the 
data and the arbitrariness in the scaling factors one can- 
not draw firm conclusions as to the significance of this. 

The very small E pp region has been studied near the 
forward direction with the COSY-ANKE spectrometer 
up to 800 MeV [H El]. The results show an energy 
dependence that suggests some influence from interme- 
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diate AY pairs. However, for kinematic reasons, at these 
higher beam energies the 7 and it peaks, which are so 
prominent in the missing-mass plot of Fig. [2j merge and 
the extraction of app — > ppj signal is much more delicate. 
Under these conditions it may be necessary to measure 
the photon in coincidence and data of this type from the 
COSY-WASA facility are currently being analyzed at 500 
and 550 MeV Hl[. 

However, it should be stressed that data in the small 
E pp region taken at well below the pion production 
threshold would also be very valuable because the un- 
certainties regarding the inclusion of the A contribution 
would then be minimized. Data with polarized beam 
and target, along with results on the differential cross 
section and analyzing power, would allow some of the 
electromagnetic multipoles to be determined explicitly 
and therefore allow one to identify defects in the models 
more clearly |33| . 

The bremsstrahlung reaction is the simplest of all the 
high-momentum-transfer reactions in proton-proton col- 
lisions at intermediate energies. Unlike cases of meson 
production, there is no need to consider the final state 
interaction of the meson with one of the protons. The 
quality of the agreement between the predictions of a 



modern bremsstrahlung model with the high-statistics 
and high-acceptance data achieved at the PROMICE- 
WASA facility is striking and should encourage further 
experimental and theoretical work in the field. 
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